Failure  Fronts  in  Brittle  Materials  and 
Their  Morphological  Instabilities 


by  Michael  A.  Grinfeld,  Scott  E.  Schoenfeld,  and  Tim  W.  Wright 


ARL-TR-3660 


September  2005 


Approved  for  public  release;  distribution  is  unlimited. 


NOTICES 

Disclaimers 

The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army  position  unless 
so  designated  by  other  authorized  documents. 

Citation  of  manufacturer’s  or  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the 
use  thereof. 


Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it  to  the  originator. 


Army  Research  Laboratory 

Aberdeen  Proving  Ground,  MD  21005-5066 


ARL-TR-3660 


September  2005 


Failure  Fronts  in  Brittle  Materials  and 
Their  Morphological  Instabilities 

Michael  A.  Grinfeld,  Scott  E.  Schoenfeld,  and  Tim  W.  Wright 
Weapons  and  Materials  Research  Directorate,  ARL 


Approved  for  public  release;  distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering 
and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to 
comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

September  2005  Final 


4.  TITLE  AND  SUBTITLE 

Failure  Fronts  in  Brittle  Materials  and  Their  Morphological  Instabilities 


3.  DATES  COVERED  (From  -  To) 

1  October  2004-30  August  2005 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Michael  A.  Grinfeld,  Scott  E.  Schoenfeld,  and  Tim  W.  Wright 


5d.  PROJECT  NUMBER 

AH84 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRD-ARL-WM-TD 
Aberdeen  Proving  Ground,  MD  21005-5066 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

ARL-TR-3660 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


14.  ABSTRACT 

This  report  outlines  the  results  of  the  effort  to  suggest  a  theoretical  method  of  describing  failure  fronts  in  glasses  and  ceramics. 
There  are  various  observations  and  experiments  showing  that  in  addition  to  standard  shock-wave  fronts,  which  propagate  with 
transonic  velocities,  other  much  slower  wave  fronts  can  propagate  within  glass  or  ceramic  substances  undergoing  intensive 
damage.  These  moving  fronts  propagate  into  intact  substance,  leaving  intensively  damaged  substance  behind  them.  They  have 
been  called  failure  waves.  In  this  report,  we  model  them  as  sharp  interfaces  separating  two  states — the  intact  and  comminuted 
states.  The  approach  is  based  on  an  analogy  between  failure  fronts  and  fronts  of  slow  combustion.  Two  main  results  are 
announced.  One  of  them  concerns  the  speed  of  a  failure  wave  driven  by  oblique  impact  of  a  brittle  target,  and  the  other 
establishes  a  criterion  for  morphological  instability  of  failure  fronts. 


15.  SUBJECT  TERMS 


ceramic  armor,  brittle  fracture,  failure  waves,  morphological  instability 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Michael  A.  Grinfeld 

a.  REPORT 

UNCLASSIFIED 

b.  ABSTRACT 

UNCLASSIFIED 

c.  THIS  PAGE 

UNCLASSIFIED 

UL 

16 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

410-306-0793 

11 


Standard  Form  298  (Rev.  8/98) 
Prescribed  by  ANSI  Std.  Z39.18 


Contents 


List  of  Figures  iv 

1.  Introduction  1 

2.  The  Main  Idea  of  the  Model  1 

3.  Formal  Statement  of  the  Problem  2 

4.  Morphological  Instability  of  the  Failure  Front  4 

5.  References  7 

Distribution  List  9 


iii 


List  of  Figures 


Figure  1 .  The  energy  balance  for  a  damaged  spring . 2 

Figure  2.  Oblique  impact  of  a  brittle  substance . 2 


1.  Introduction 


Extensive  experiments  with  glasses  and  brittle  ceramic  materials  made  by  different  groups  in 
Britain,  Russia,  and  the  United  States  in  the  1990s  (1-8)  proved  that  failure  waves  exist.  This 
conclusion  is  backed  by  other  observations  in  geomechanics  and  engineering  (9-12). 

The  problem  of  failure  waves  demands  significant  progress  of  the  relevant  experiment  (theory 
and  numerical  modeling).  These  three  pillars  grow  simultaneously  and  in  close  interaction  with 
each  other.  The  experimenters  find  an  analogy  of  such  waves  with  solid-solid  phase 
transformation  fronts  (13)  or  with  fronts  of  slow  combustion  (14).  These  two  analogies  are  not 
antagonistic.  Modeling  of  both  phenomena  includes  changes  in  the  thermodynamic  potentials  of 
the  solid  states  involved.  Although  damage  is  definitely  not  a  solid/solid  phase  transformation, 
both  problems  have  the  same  roots — it  is  the  minimization  of  accumulated  energy  by  big 
changes  in  the  microstructure.  For  us,  the  analogy  with  combustion  looks  somewhat  more 
relevant  and  appealing. 

Our  model  is  based  on  the  analogy  with  the  simplified  theory  of  slow  combustion.  In  the  theory 
of  slow  wave  with  phase  transfonnations,  such  an  approach  was  quite  successfully  applied  about 
two  decades  ago  in  geophysics  and  celestial  physics  (15)  and  in  low-temperature  physics  (16). 
Certain  difficulties  of  the  simplified  theory  of  phase  transformation  waves  have  been  discussed 
in  Grinfeld  (17)  and  Glimin  (18).  The  discussion  given  there  is  equally  relevant  for  modeling  of 
failure  waves.  According  to  the  simplified  theory,  the  internal  wave  “structure”  of  the  wave 
front  is  ignored,  and  the  failure  front  is  treated  as  a  mathematical  surface.  (Relevant  discussions 
can  be  found  in  the  classical  monographs  in  Landau  and  Lifshitz  [79]  and  Courant  and  Friedrichs 
[20].) 


2.  The  Main  Idea  of  the  Model 


To  justify  the  main  idea  of  the  energetic  approach,  let  us  consider  an  elementary  spring-like 
system  like  that  shown  in  figure  1 .  The  right  edge  of  the  rod  is  fixed.  The  left  edge  is  under  the 
action  of  an  axial  loadT3  .  The  load  is  large  enough  to  cause  a  considerable  damage  in  the  rod. 
The  rod  is  divided  by  the  failure  front  F  into  the  intact  domain,  shown  in  white,  and  the 
damaged  domain,  shown  in  black.  For  the  sake  of  simplicity,  we  assume  that  the  intact 
substance  has  the  Young’s  modulus  A, ,  whereas  the  weaker  damaged  substance  has  the  Young’s 
modulus  Ed  <  Ei .  The  total  energy  of  the  model  system  is  assumed  to  be  equal  to 

EM=AP+(Eyj/2+eb)x+E,e;(L-x)/2,  (1) 
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Figure  1.  The  energy  balance  for  a  damaged  spring. 

where  A  is  the  displacement  of  the  left  end  of  the  rod,  si  and  ed  are  the  defonnations  in  the  two 
domains  of  the  combined  rod,  and  eb  is  the  energy  of  the  broken  bonds  per  unit  length.  In  the 
quasistatic  evolution,  we  obviously  get  the  two  following  equations: 

A  =  etx  +  ed{L-x),  (2) 


and 


p  =  Eded 


=  Eiei. 


(3) 


Using  the  last  two  equations,  we  can  rewrite  the  total  energy  of  the  system  as  follows: 


J  total 


et-(El-E7')P2/ 2 


x  +  P2L/2E;. 


(4) 


The  last  fonnula  implies  that  at  a  sufficiently  large  load,  P ,  such  that  (E7-E;')P2>leb,  tlie 
growth  of  the  length  x  of  the  damaged  zone  becomes  energetically  favorable. 


3.  Formal  Statement  of  the  Problem 


We  limit  our  study  with  two-dimensional  propagation  and  consider  an  initially  resting  unifonn 
half  plane  x  >  0  experiencing  an  impact  at  x  =  0  by  the  oblique  force  P  (figure  2). 
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Figure  2.  Oblique  impact  of  a  brittle  substance. 
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Within  each  of  the  domains,  the  energy  densities  eint  and  edam  per  unit  mass  of  the  intact  and 
damaged  states,  respectively,  are  given  by  the  following  formulas: 

Peint  0 Um,n  )  =  0 1  /2)4f  >  (5) 

and 

Pedam  (Um,n  )  =  ( l/2)CdlU,,jUkJ  +  &  ’  (6) 

where  umn  is  the  displacement  gradients,  c^f  and  cf^m  are  the  elasticity  tensors  of  the  two  states, 
and  p  is  the  original  mass  density  (a(.  is  the  symbol  of  differentiation  with  respect  to  the  spatial 
coordinates  x  ).  Limiting  ourselves  with  the  approximation  of  linear  elasticity,  we  ignore  all 
effects  of  mass  density  change.  The  positive  constant  qh  takes  into  account  the  energy  required 
to  produce  various  defects  (interfaces,  vacancies,  shear  bands,  holes,  etc.)  distributed  within  unit 
mass  of  the  bulk  of  damaged  substance.  This  term  is  analogous  to  the  constant  used  in  the  theory 
of  slow  combustion,  which  takes  into  account  the  energy  release/consumption  due  to  chemical 
reactions.  (See  Grinfeld  and  Wright  [21  ]  for  a  detailed  discussion  of  the  model.) 

In  addition  to  appropriate  initial  and  (external)  boundary  conditions,  the  master  system  includes 
the  bulk  equilibrium  equation  within  each  of  the  bulk  domains  (equation  7),  the  displacement 
continuity  equation  across  the  failure  front  (equation  8),  and  the  traction  continuity  across  the 
failure  front  (equation  9). 
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where  p"  =  pde(umn)/dujj  is  the  stress  tensor  and  n  is  the  unit  normal  to  the  failure  front. 

The  last  equation  (equation  10)  across  the  failure  front  describes  the  kinetics  of  failure  as 
follows: 


c  =  -K[pJ-k']_n.nk ,  (10) 

where  c  is  the  velocity  of  the  failure  front,  the  tensor  p'k  is  defined  as  p'k  =  edJk  -  p  '  p"  (dik  +  uik )  — 
it  plays  the  same  role  as  the  scalar  Gibbs  chemical  potential  p  of  a  liquid  substance — and  K  is  a 
positive  (kinetic)  constant  or  function  with  dimension  [velocity]-1.  We  assume  explicitly  that  the 
displacements  and  traction  remain  continuous  across  the  interface.  This  is  not  the  only 
reasonable  option.  Another  reasonable  option,  especially  when  dealing  with  pulverized  states, 
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would  be  the  model  of  a  friction-free  interface  with  discontinuous  displacements.  In  this  case, 
the  last  kinetic,  constitutive  equation  (equation  10)  should  be  modified  as  well  (a  similar  system 
was  analyzed  in  Parshin  [16]  and  Grinfeld  [77]  in  the  context  of  phase  transformations). 

The  system  (equations  7-10)  allows  piecewise  linear  solutions  of  the  following  forms: 

u\(x,z,t)  =  d'+t  +  a'+z  at  x>  ct  ,  (11) 


and 


u‘_ (x,  z,  t)  =  d'_t  +  a1  z  at  x  <  ct  , 


(12) 


where  d'±,  a\ ,  and  c  are  certain  constants.  The  sign  ±  marks  the  quantities  related  to  the  intact 
(damaged)  state.  The  constant  c  gives  the  speed  of  the  failure  front  across  which  the 
displacement  gradients  a\  suffer  the  finite  discontinuities.  The  solution  just  mentioned  allows 

considering  the  quasistatic  problem  oblique  impact  (loading)  when  the  applied  force  has  the 
following  components: 

P±=n,  (13) 


and 


Pl=T  . 


(14) 


In  this  case,  the  system  (equations  7-10)  leads  to  the  following  formula  of  the  velocity  of  the 
failure  front: 


2 

P  _K 
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1 

1  > 
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f---l 
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K  2 
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2 

V  d-  d+  y 

(15) 


4.  Morphological  Instability  of  the  Failure  Front 


In  order  to  explore  morphological  stability  of  the  piecewise  linear  solution,  we  present  the 
elastic  displacement  u’±(x,z,t)  and  the  speed  of  the  interface  c(z,t)  in  the  following  forms: 

u± (x, z,t)  =  u°l (x, t)  +  u±(x,z,t) ,  (16) 


and 


c(z,t)  =  c°  +c  (z,t) ,  (17) 

where  ujx,z,t)  and  c(z,t)  are  small  disturbances.  These  equations  should  be  substituted  in  the 
bulk  equations  and  in  the  boundary  conditions  (equations  7-10),  which  should  then  be  linearized  with 
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respect  to  the  small  disturbances.  We  then  look  for  the  solutions  of  the  linearized  system  in  the 
following  forms  (k  is  the  in-plane  wave-number  and  //  is  the  rate  of  growth): 

u± ( x , z,t)  =  W±(x  —  c°t'J ekz+v ,  (18) 

and 

£{z,t)  =  Seik:+nt ,  (19) 

where  the  functions  Wl  ( X )  should  decay  exponentially  at  X  — »  ±  oo  .  The  dispersion  equation 
for  q  is  too  lengthy  to  be  presented  here.  It  becomes  relatively  compact  and  instructive  in  the 
limit  of  an  incompressible  damaged  state  (v_  =  1/2  )  and  v+  =  1/3  (v±  =  XJ2{X.  +  74))  when  it 
reads  as  follows: 


p  }  T2  M+(m+~M-)2  4p  +  3/4 

8 Kk  5 pi  +87474  +3pl  11“  p  5 pi  +8p  p++3pi 

The  wave  front  is  morphologically  unstable  if  the  linearized  master  equation  has  solutions  of  the 
previous  form  and  with  the  rate  q  corresponding  to  exponential  growth  in  time.  The  last  formula 
shows  that  shear  stresses  T  play  a  stabilizing  role  (similar  to  the  case  of  morphological 
instabilities  of  solid-solid  phase  interfaces  [77]).  In  addition,  at  T  =  0,  the  failure  front  is 
morphologically  unstable  in  the  most  interesting  case  when  the  shear  modulus  p  of  the  damaged 
state  is  less  than  the  shear  modulus  p+  of  the  intact  state.  Speaking  very  schematically,  the 

instability  has  a  simple,  physical  meaning.  It  means  that  penetration  of  fingers  of  damaged 
material  into  intact  material  is  the  fastest  way  of  releasing  accumulated  elastic  energy  from  the 
system. 

The  morphological  failure-front  instability  has  the  potential  to  explain  appearance  of  fingering 
cracks  on  disintegration  fronts  in  Prince  Rupert  drops  reliably  documented  in  experiments  of 
Chandrasekar  and  Chaudhri  (4).  This  interpretation,  however,  requires  further  (numerical) 
studies  of  deeply  nonlinear  stage  of  the  instability.  Another  appealing  possibility  of  the 
instability  applications  is  to  the  failure-front  waviness  that  was  detected  in  the  numerical 
modeling  of  Resnyansky  et  al.  (22).  These  authors  use  a  formally  different,  but  conceptually 
close,  model.  Appearance  of  the  waviness  in  their  simulations  could  be  not  a  numerical  artifact 
but  a  manifestation  of  the  instability  discussed  here. 

In  conclusion,  we  have  presented  a  simple  model  for  a  failure  front  based  on  an  analogy  with 
slow  combustion  and  a  fonnula  for  the  failure  wave  generated  by  oblique  impact  on  brittle 
material.  We  then  demonstrated  that  under  rather  general  assumptions,  the  flat  failure  front  is 
morphologically  unstable.  The  last  conclusion  shows  an  additional  analogy  between  failure 
fronts  and  slow  combustion  fronts  (but,  of  course,  the  mechanisms  of  destabilization  are  totally 
different).  The  morphological  instability  discussed  here  may  shed  some  light  on  the  appearance 
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of  corrugations  on  failure  fronts  and  crack  bifurcations  as  observed  in  experiments  on  static 
indentation  of  brittle  materials,  the  appearance  of  radial  cracks  in  dynamic  experiments  with 
penetration  of  projectiles  through  brittle  materials,  the  appearance  of  corrugations  on  the 
disintegration  fronts  in  Prince  Ruppert  drops  ( 4 ),  and  the  appearance  of  failure  front  waviness 
detected  in  numerical  simulations  (22). 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRD  ARL  RO  EN 
A  RAJENDRAN 

RESEARCH  TRIANGLE  PARK  (ARO) 
DURHAM  NC  27703 

13  JOHNS  HOPKINS  UNIVERSITY 
DEPT  OF  MECHL  ENGRG 
F  ZHOU 
K  VOLOKH 
T  HUFNAGEL 
Z  H  AIT  AO 
EMA 

J-F  MOLINARI 
K  J  HEMKER 
K  TRAMESH 
Q  WEI 
B  PALIWAL 
G  WENG 
S  C  DANFORTH 
W  TIM 

LATROBE  122 
3400  N  CHARLES  ST 
BALTIMORE  MD  21218 

ABERDEEN  PROVING  GROUND 


51  DIRUSARL 

AMSRD  ARL  Cl  HC 
D  GROVE 
AMSRD  ARL  WM 
M  GREENFIELD 
T  WRIGHT 
AMSRD  ARL  WM  M 
J  BEATTY 
J  MCCAULEY 
AMSRD  ARL  WM  MB 
KCHO 
R  DOWDING 
L  KECSKES 
M  ST  AKER 
J  SWAB 
R  WOODMAN 
AMSRD  ARL  WM  MC 
R  ADLER 
J  LASALVIA 
AMSRD  ARL  WM  MD 
J  ADAMS 
J  CAMPBELL 
B  CHEESEMAN 
ECHIN 
K  DOHERTY 
C  FOUNTZOULAS 


G  GAZONAS 
S  GHIORSE 
G  GILDE 
P  HUANG 
P  PATEL 
WROY 
T  SANO 
C  YEN 

AMSRD  ARL  WM  TA 
M  BURKINS 
W  GOOCH 
V  HERNANDEZ 
C  HOPPEL 
E  HORWATH 
D  KLEPONIS 
C  KRAUTHAUSER 
B  LEAVY 
D  MACKENZIE 

AMSRD  ARL  WM  TC 
R  COATES 
M  FERMEN  COKER 
L  MAGNESS 
B  SCHUSTER 
C  WILLIAMS 

AMSRD  ARL  WM  TD 
D  CASEM 
D  DANDEKAR 
K  IYER 
H  MEYER 
M  RAFTENBERG 
E  RAPACKI 
M  SCHEIDLER 
S  SCHOENFELD 
S  SEGLETES 
T  WEERASOORIY A 
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